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Abstract 

We perform a detailed comparison of long range rapidity correlations in the Color Glass Con- 
densate framework to high multiplicity di-hadron data in proton-proton and proton-lead collisions 
from the CMS, ALICE and ATLAS experiments at the LHC. The overall good agreement thus far 
of the non-trivial systematics of theory with data is strongly suggestive of gluon saturation and the 
presence of subtle quantum interference effects between rapidity separated gluons. In particular, 
the yield of pairs collimated in their relative azimuthal angle A(f) ~ 0, is sensitive to the shape of 
unintegrated gluon distributions in the hadrons that are renormalization group evolved in rapidity 
from the beam rapidities to those of the measured hadrons. We present estimates for the collimated 
di-hadron yield expected in central deuteron-gold collisions at RHIC. 



1 



I. INTRODUCTION 



In two recent papers [J, 0], we argued that data on two particle correlations in high 
multiplicity proton-proton and proton-nucleus collisions from the CMS collaboration [3|, 
3| provided strong evidence for gluon saturation and the Color Glass Condensate (CGC) 
effective field theory (EFT) jsj describing this phenomenon. These correlations, which are 
long range in the relative rapidity Arj between pairs of charged hadrons, show an unexpected 
"nearside" coUimation A(/) ^ in their relative azimuthal angle. This unusual coUimation 
is called the "ridge" due to its structure in the Ar; — plane. Reviews of this ridge effect 
can be found in (g], 0] ■ 

In the CGC effective theory, the nearside coUimation is obtained from QCD graphs called 
"Glasma graphs"; for high occupancy gluons (with transverse momenta kj_ < Qs, the sat- 
uration scale), Glasma graphs are enhanced by a^^, a factor of ~ 10^ for typical values of 
the probed QCD fine structure constant as- In the power counting of the EFT, the effect of 
gluon saturation on Glasma graphs ensures they provide a significant additional contribution 
in high multiplicity events to "di-jet" QCD graphs. The latter are kinematically constrained 
to provide an "awayside" back-to-back coUimation peaked at A(j) ^ n but do not provide a 
significant nearside coUimation. 

The importance of Glasma graphs was first discussed in jsf and the formalism developed 
in [qI, [io[. It was first postulated as an explanation of the high multiplicity CMS proton- 



proton ridge in and a quantitative description of the nearside coUimated yield obtained 
in (H. 

In the first of a current series of papers pL|] , the description of long range di-hadron corre- 
lations in high multiplicity events was significantly developed by considering both nearside 
and awayside collimated contributions. In the former case, the Glasma graphs provide the 
dominant contribution, while the awayside receives contributions from both Glasma and 
back-to-back QCD graphs. The latter, in the high energy kinematics of the LHC experi- 



ments, is described in the CGC EFT by BFKL dynamics [13|, [IJ]. We showed in this study 
that the BFKL dynamics, which generates gluon emissions between gluons that fragment 
into triggered hadrons, does well in describing the awayside spectra. The description is 



significantly better than PYTHIA-8 [3[, and 2 — )• 4 QCD graphs in the Quasi-Multi-Regge- 



Kinematics (QMRK) [15|, ll6[, both of which overestimate the awayside yield, especially at 
larger momenta. 

In the second paper in this series j2|, we applied the Glasma -l-BFKL CGC framework 
to describe first p+Pb data on the ridge obtained at ^/s = 5.02 TeV/nucleon by the CMS 
collaboration ^,]. The CMS p+Pb data had the following striking systematic features: i) a 
strong dependence of the ridge yield on the number of charged particle tracks iVtrack, with a 
significantly larger signal than in p+p for the same A^track, ii) a stronger pt dependence than 
in p+p for the same large A'track, and iii) a nearside coUimation for large iVt^ack comparable 
to the awayside for the lower px = p%^^ = p^^ di-hadron windows. In j2|, we showed that 
all these systematic features could be explained in the CGC framework as a consequence of 
a remarkable quantum intereference effect in the production of correlated gluons. 

Subsequently, both the ALICE [13] and ATLAS [3] collaborations have presented their 
di-hadron correlation results from the first LHC p-|-Pb run. The ALICE experiment has 
an acceptance in A?] of |A?7| < 1.8, while the ATLAS experiment has an acceptance of 
2 < lAr/l < 5, close to the CMS acceptance of 2 < |A?7| < 4. In addition to the LHC 
results, the PHENIX collaboration at RHIC have reanalyzed their deuteron-gold data at 
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200 GeV/nucleon and have extracted a ridge signal in very central events [19[. All three of 
the experiments show that when the two particle yield in peripheral collisions is subtracted 
from the central events, a dipole structure remains that is long range in rapidity. This is 
precisely what one would anticipate in our Glasma+BFKL graph scenario, because the latter 
has a weak dependence on centrality, and the former has shape that is symmetric around 
A0 = 7r/2. In this work, we will address these recent analyses, and show that they can be 
reproduced in the CGC framework with a common set of parameters. 

This work is organized as follows. In Section 2, we will very briefly introduce our frame- 
work. This discussion will hew closely to those in Refs. with many details to be found 
in those papers and the papers cited therein. In section 3, we will first discuss results for 
proton-proton collisions reported by CMS. In addition to the "matrix" of the collimated 
associated yield per trigger as a function of for varying windows in p)^'^, p'^'^ for the 
highest multiplicity events, we will show results for different multiplicity bins as well. We 
will comment on some open issues here that may impact the interpretation of the p-|-Pb 
data. We next show a comparison of our results for p+Pb collisions for the associated yield 
per trigger to data from CMS, ATLAS and ALICE. We also present estimates for the associ- 
ated yield in deuteron-gold collisions at RHIC, which is likely to be obtained soon following 
preliminary results presented by the PHENIX collaboration [l^. In the final section, we 
summarize our results and discuss their implications, as well as comment on alternative 
interpretations of the data. An appendix outlines the different normalization procedures 
followed by the different experiments, and how these can be related to our analysis and to 
each other. 



II. REVIEW OF LONG RANGE RAPIDITY CORRELATIONS IN THE CGC EFT 

In this section, we will review the expressions discussed in our previous papers. Except 
for a few clarifying details, the discussion is very similar to that in [ll, (3]. In Fig. ([1]), as in 
the previous two papers, we provide a schematic sketch of the Glasma and BFKL graphs. 
Both of these are connected QCD graphs, and their relative power counting is given by the 
CGC EFT. As shown in the figure, the Glasma graphs give a collimated contribution in Acf) 
that is mirror symmetric about Acf) = 7r/2. The "di-jet" contribution is peaked back-to-back 
around A0 = tt, and gives a negligible contribution on the nearside at A0 = 0. 

The collimated contributions from all the Glasma graphs can be compactly written as [2I, 
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(1) 



where we define 

Di = '^l.iVp, kT)<^A2{yp, Pt - kr) ["^A^iVg, Qt + kr) + $^2(2/9, Qt - ky)] , 

D2 = iVq, kT)$Ai iVp, Pt - kr) {Vg, qr + ^t) + {Vg, qr - k^)] . (2) 

These four terms, called the "single diffractive" and "interference" graphs in [^], constitute 
the leading pr/Qs behavior. Also included is the next order correction in pr/Qs where we 
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FIG. 1. Anatomy of di-hadron correlations. The glasma graph on the left illustrates its its schematic 
contribution to the double inclusive cross-section (dashed orange curve). On the right is the back- 
to-back graph and the shape of its yield (dashed blue curve). The grey blobs denote emissions all 
the way from beam rapidities to those of the triggered gluons. The solid black curve represents 
the sum of contributions from glasma and back-to-back graphs. The shaded region represents the 
Associated Yield (AY) calculated using the zero-yield-at-minimum (ZYAM) procedure. Figure 
from ref. till. 



have|20[ Ai = 5^{pt + qr) [If + + 2X|], such that 

(ki_L ■ Pt - ^l±f 



^1= $Ai(yp,kix)$A2(2/g,Pr-kix) 

^kix kf^ [pt - kix) 



2 ' 



^2= $Ai(yp, ki_L)$A2(2/g, Pt - k 

^3 = / $Ai(yp,kix)$A2(?/g,PT - k 



|ki_L X Pt 



2 



k^^ (Pt -ki^'f ' 



(ki_L ■ Pt - ki±) |ki_L X pt\ 
K± (Pr -ki±? 

The other contribution, A2, in Eq. ([1]) can be expressed as 

^2= / '^'Ai(yp,ki_L)$Ai(?/p,k2_L)$A2(2/q,PT - kix)$A2(yq,qT + ki_L) 

(kix ■ Pt - kl±) (k2± ■ Pt - \4±) + (ki± x pt) (ksx x pt) 



kL (Pt - k 



l-L 



^ (kix • Qt - kf±) (k2± ■ Qt - k^±) + (ki± x qt) (kax x qt) 

kL (qT + ki±f 

where k2_L = pt — Qt — ki^. The above expressions are the result of including all combina- 
torial combinations of graphs represented by the Feynman diagram to the left in Fig. [1] The 
combinatorics is a result of different ways of averaging over strong color sources between the 
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amplitude and complex conjugate amplitude in both projectile and target. It is important 
to note that Eq. ([1]) represents genuine quantum interference contributions, the structure of 
which were first outlined in [8|. 

In equations[T]through[3]the only function (besides the one loop running coupling constant 
as) is the unintegrated gluon distribution (UGD) per unit transverse area 

<^A{y,k^) = ^^^^ r dr^r^Uk^r^)[l-TA{y,r^)f (4) 
^as Jo 

where Ta is the forward scattering amplitude of a quark- ant iquark dipole of transverse size r± 
on the target A; it, or equivalently, the UGD, is a univers al q uantity that can be determined 



by solving the Balitsky-Kovchegov (BK) equation |2l|, |22| as a function of the rapidity 
y = log (xq/x). Specifically, what we use for the UGDs is the rcBK equation, which includes 
all leading logs in x (LLx) contributions to the UGDs + running coupling next-to-leading- 
logs (NLLx) effects via the Balitsky prescription [23|. The forward scattering amplitude 
TA^y, r±) at the initial scale x = xq is a dimensionless function of rj_Ql, where Qo is a non- 
perturbative scale at the initial rapidity. The saturation scale Qs, defined as the transverse 
momentum defining the peak value of $ on the l.h.s of eq. (jl]), is typically a larger scale 
even at the initial rapidity, and grows rapidly via the rcBK renormalization group equation 
with rapidity. In the rcBK equation, different impact parameters in the proton/nuclear 
target are modeled by varying Qo- The minimum-bias (median impact parameter) value we 
choose for the proton Ql = 0.168 GeV^ (corresponding to a Qs ~ 0.7 GeV in the adjoint 
representation at the initial rapidity), is the value that gives a best fit to deeply inelastic 
electron-proton scattering data from HERA jl^ . 

In addition to Qq, Eq. ([1]) has three free parameters; the first being the transverse over- 
lap area S± is fixed separately for p-|-p, d+Au and p+Pb collisions and will be discussed 



in detail later. Once 5*^ is fixed [25| for a given system, the centrality dependence is con- 
trolled entirely though the choice of initial saturation scales for the projectile and target. A 
second parameter, held fixed to the same value in both p+p and p+A collisions, is the non- 
perturbative constant ( = 1/6 specifying the correction to the factorized UGD description 
due to soft multigluon interactions. It is independently constrained by fits to empirical p+p 
multiplicity distributions within the k± factorization approximation for multi-gluon produc- 
tion [26|, |27| and by real time classical Yang-Mills computations We note further 



that the collimated structures seen in the perturbative classical computations persist in the 
full non-perturbative classical results, thereby lending confidence that the latter primar- 
ily renormalize the amplitude of the former. The third parameter -ft'ciasma we will discuss 
shortly. 

The framework of Glasma graphs is based on the factorization theorems for "dense- 
dense" systems [sot] derived in [9'], which include leading log corrections to all orders in 
perturbation theory (so called LLx approximation) as well as all leading multiple scattering 
contributions l3l|. As the full expression is very cumbersome, a Gaussian truncation is 
employed in [lO|, where Eq. ([T]) was first derived. The Gaussian truncation was shown in 
to be a very good approximation to the full JIMWLK evolution. In addition, it is 



assumed that Qs < ky, in order to obtain the expression in terms of UGDs (unintegrated 
gluon distributions). We emphasize that the resulting expression cannot be interpreted 
simply as the product of UGDs with matrix elements, but combines LLx contributions to 
each. 

We now consider the double inclusive distribution from the back-to-back BFKL graphs 
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33| shown in Fig. [H The double inclusive multiplicity can be expressed as [l5|, |34 



K^m (5) 



X 



/ / $A(a;i,ko±)$B(x2,k3_L) ^(ko± - pT,k3_L + qT,2/p - 



where Q is the BFKL Green's function 



11 — /'+°° / \ 



is the BFKL eigen- 
le Gamma function. 



Here Cp = {N^ - l)/2N,, uj{iy,n) = -2a, Re ^ + - ^(1) 

value, where "^{z) = dlnr{z)/dz is the logarithmic derivative of t 

Further, we have a, = NcOs (\/Qa±Qb±) /tt and 6 = arccos ( r^Hr i 

In the description of the away-side jet in the above BFKL framework, the UGD evolution, 
as for the Glasma graphs is described by the rcBK equation. NLLx corrections to the back- 
to-back graphs have been computed in jl^]. It was demonstrated there that the NLLx 
correction to the A0 independent pedestal is a large one (a factor 2 to 3). However, the 
NLLx contribution to the collimated (cos(A0)) and (cos(2A(/))) moments (as also confirmed 
in jlsf), which are the quantities of interest here, is 10-30%; further, we expect our inclusion 
of running coupling in Eq. ([6]) will potentially account for a good fraction of this correction. 
Based on the results in these works, it is reasonable to conclude that the BFKL contribution 
to the collimated yield has 10-30% uncertainties. 

As shown in Fig. ([1]), Eq. (|5]) gives a collimated A0 contribution exclusively on the away 
side, peaked at A0 = vr, while Eq. ([1]) gives a "dipole" cos(2A0)-like contributions with 
maxima at and vr. It should be noted that the behavior of both sets of graphs for any 
given window in and rapidity is strongly influenced by the evolution of the UGDs which 
are common to both and provide a non-trivial constraint on the relative contributions of each 
to the description of data. It's the interplay between these contributions with varying Qo 
in projectile and target that describes the systematics of the proton-proton and proton-lead 
data, that we shall discuss quantitatively shortly. 

The K-factors used throughout this work, -f^Giasma and -R'bfkl, are introduced to take into 
account not only uncertainties in higher order computations but also acceptance corrections 
and uncertainties in the choice of fragmentation functions. For simplicity, we will take them 
to be equal to each other in each process, but will use different values, as stated at the ap- 
propriate juncture, for p+p and p-|-A. There is no reason a priori why all these should be the 
same in p-|-p and p-|-A. Also, modulo a better understanding of the multiplicity distribution 
in p+A, some of the uncertainties in K factors could be absorbed in the initial saturation 
scale Qo or vice versa, corresponding to slightly different number of participants. In addi- 
tion to data on multiplicity distributions in the rapidity window of interest, uncertainties 
on fragmentation functions in particular can be constrained by forthcoming data from the 
LHC on single particle spectra at forward rapidities. 
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FIG. 2. The nearside yield per trigger as a function of (specific to the acceptance of the CMS 

experiment) for 1 < < 2, for = "p^^ = p'^'^ . The open circles (brown) are the computed 
Glasma graph yield for increasingly rare proton initial saturation scales in multiples -A^part°" of 
Qo = 0.168 GeV^. The filled circles (blue) are proton-proton data at y/s = 7 TeV from the CMS 
collaboration. Each of the p+Pb curves corresponds, for fixed QolP^o^on), to the increasing yield 
with a larger number of participants in the nucleus (in increments of two, from N^^j.^ = 6 to 
-^part — 22), where N^^^.^ denotes the initial saturation scale in the Pb nucleus through the relation 
Q§(lead) = N^^^^ ■ 0.168 GeV^. The three curves with open symbols are new computations for 
higher multiplicities. In this case the number of participants in the nucleus increases in increments 
of two up to A^part = 26. The red squares are CMS data for proton-lead collisions at y/s = 5.02 
GeV/nucleon. Similar curves are obtained for the -/V^^f™*^ appropriate in the acceptance of the 
other experiments. 



III. QUANTITATIVE STUDY OF DI-HADRON CORRELATIONS AT WIDE RA- 
PIDITY SEPARATIONS 

For the analysis of the LHC proton-proton and proton-lead data, as well as the RHIC 
deuteron-gold data, we must be able to make a reasonable estimate of the centrality class 
based on the total charge particle multiplicity. In the CGC framework, the single inclusive 
gluon distribution is defined as [12i] 

dA^i asNc S± f . 

dP^ = A^mh^Wr JJ^^(y^^^T)^AM,PT-kT), (7) 
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FIG. 3. Long-range (2 < \Ar]\ < 4) per-trigger-yields [l / Ntngd? N / dA(f^ of charged hadrons as a 
function of \A(j)\, for p+p collisions at ^/s = 7 TeV. Data are from the CMS collaboration. The lower 
(upper) curves indistinguishable in some windows, correspond to the following: i) N°^^^'^ < 35: 
N^Tr = 1,2 ii) 35 < iV°f- < 90: = 3,4 iii) 90 < iV°f- < 110: N^lT = 4,5 iv) 

110 < iV°f < 150: iVPart"" = 5,6 v) iV°f > 150: N^^"'^ = 7,8. 

where the $'s are the UGDs defined previously. The number of charged hadron tracks is 
then defined as 

/+?/acccpt-J/shift r i^AT 
dV d'vT-r-^{VT). (8) 

where ?/shift = 0.465 is the shift in rapidity in the center-of-mass frame in asymmetrical p+Pb 
collisions towards the lead fragmentation region, j/accept is the maximal laboratory rapidity of 
a given detector and pT,min is the minimal transverse momentum of measured charged tracks 
in the detector [sg!] . Since Eq. ([7]) corresponds to the single inclusive gluon multiplicity, and 
since the bulk of particle production is soft, we have to introduce a gluon liberation factor 
Hg, which specifies, assuming parton-hadron duality, the conversion of gluons to charged 
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FIG. 4. Long-range (2 < [Ar^l < 4) per-trigger-yields [l/Ntj-igd?N/dA(fj of charged hadrons as a 
function of \ from high multiplicity (-/V^°k'™'^ — 1^^) PP collisions at ^/s = 7 TeV. Data are from 
the CMS collaboration. The lower (upper) curves, indistinguishable in some windows, correspond 
5,6. 



to 



hadrons. In Eq. ([H]), the combination of the transverse overlap area S± times Kg is separately 
fixed for p+p and p+Pb to give the best description of the data. For minimum bias proton- 
proton collisions {Qq = 0.168 GeV^ in both protons) this corresponds to A^°rk™^ ~ l^. 

This value, {KgS±), is subsequently held fixed to determine N°^™^ as Ql (the saturation 
scale at the initial rapidity for x evolution) in both the proton and lead nucleus is varied. 
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Admittedly, the constant {KgS±) provides a lot of freedom in the centrality selection, but it 
can only be constrained as further data on single inclusive quantities in the rapidity ranges 
of interest become available. The uncertainty in the value of {KgS±) results in a rescaling 
of the X-axis of Fig. |21 Neither Kg or S± enter into the computation of the associated yield. 
As discussed later, the interaction cross section is on the order of the size of the proton and 
not the nucleus. 

A part of the analysis of the di-hadron data requires the calculation of the number of 
trigger particles, defined here as 

^tng = / dJ d'pT dz y . (f ) , (9) 

where PT,min) PT,max. dcuotc the width of the pt window wherein triggered particles are 
selected. Since the triggered hadrons are semi-hard, one takes account of the possibility that 
they were generated by the fragmentation of higher pt gluons by fragmentation functions. 
These are chosen, as in to be the NLO KPP parametrization js^ of the fragmentation 
function of gluons to charged hadrons. 

The double inclusive multiplicity of charged hadrons is computed as 

d? N /"J/acccpt— ?/shift 



dA^ 



/2/acccpt— ?/shift 
drjp drjq A {rjp, rjg) (10) 



^accept yshift 

dp^ 



20 



Z2 d'^pTd'^qrdripdriq \zi' Z2 



Bounds on the range of the trigger and associated hadron momenta are denoted respectively 
as p™™'^™'^'') and Likewise, Ari^iy^{Arimax) denote the pseudo-rapidity separation 

between the measured hadrons for a given detector j38|. The acceptance 
A (rjp, rig) takes into account the acceptance of the uncorrelated background. The different 
treatments of this acceptance function between the ATLAS, CMS and ALICE experiments 
are discussed in the Appendix. The PHENIX analysis is identical to that of ATLAS. 

The coUimated associated yield is computed using the Zero-Yield-at-Minimum (ZYAM) 
procedure, 

... , 1 1 /"^"^"'"r .Jd'^N d^N 
Assoc. Yield = — — / dA(f) 



A^trig Jo \ dA(p dA(p 



(11) 



A0„ 



where A^min. is the angle at which the two particle correlation strength is minimal. An 
important point to note is that the transverse overlap area S±_ cancels out between the 
numerator and denominator in the r.h.s eliminating a source of uncertainty in di-hadron 
spectra. 

After these preliminaries, we are now ready to discuss our results. In Fig. ([2]), we plot the 
integrated associated nearside yield per trigger (obtained from Eqs. (fTOj) and ( [TT]) ) versus 
N°^™^ as determined in Eq. (|H]) for 1 < < 2, for px = p%^^ = p^^- The associated yield 
in p-|-p collisions in the computation is shown by the open brown circles that are connected 
by dashed brown lines. These open circles correspond to results in integer multiples of 
Ql = 0.168 GeV^, which as noted previously is the saturation scale at the initial xq in fits 
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of the rcBK equation to inclusive deeply inelastic scattering data at HERA. The blue filled 
circles correspond to the CMS proton-proton data for the nearside associated yield in the 
different iV°r^™'^ windows specified by the collaboration. We observe that they lie nicely on 
this curve. This then helps us identify the range in Ql (proton) that matches the N^^'""^ 
centrality selection in the experiment to be discussed shortly in comparisons to the detailed 
matrix of the collimated yield versus A0. 

The sole inputs for the trajectories shown in Fig. (j2]) are Qq (proton) =A^p™t°'^ ■ 0.168 
GeV^ and (5g(lead) =N^^^^ ■ 0.168 GeV^. Each of the curves corresponds to a fixed Ql in the 
proton of 0.168 — 0.672 GeV^ (or Npl°l°^=l-A) representing estimates of these quantities from 
median ( "min. bias" ) impact parameters in the proton to the very central impact parameters 
respectively that are triggered in high multiplicity events. The trajectories corresponding to 
each of these proton Qq show how the yield increases with a larger number of participants 
in the nucleus. Because the CMS p+Pb data have the same Arj and centrality selections, 
we plot these as well. The message one draws from interpreting this figure is that one is 
not only accessing rarer Fock configurations (at a given impact parameter) in the nucleus 
with increasing N^^^^ but also rare Fock states in the proton represented by the increasing 

Fig. ([2]) captures the essence of variations in the Glasma yield with A^t"^™*^ and iV^Jj.j.. It 
shows clearly that the yield in central p/d+A collisions is significantly enhanced relative to 
the yield in p+p collisions for the same A^"^™*^. The underlying physics behind these curves 
is the quantum interference of the UGDs and the sensitivity to the spectrum of gluons in 
the projectiles. A detailed analysis of these systematics was performed in j2| and we refer 



the interested reader to the discussion there 39 



We will now move forthwith to a comparison to data on the collimated nearside and 
awayside yields in proton-proton and proton- nucleus collisions at the LHC. In [l|, we pre- 
sented a comparison to the high multiplicity CMS data for 7 TeV proton-proton collisions. 
Subsequently, the CMS collaboration presented a detailed matrix in several centrality win- 
dows of the associated yield as a function of Acj) jl] . We also realized that the normalization 
of the CMS acceptance was different from what we assumed it to be-see the appendix for a 
detailed discussion. 

With this additional information, we have performed here a reanalysis of the CMS proton- 
proton data. The results are shown in Figs. ([3]) and (jlj) and correspond to pt = Pr^ = Pt'^- 
The comparison to data in Fig. is shown for the first time while a comparison to (jl]) 
was shown previously. The key difference to our previous comparison is that our fits for 
p+p are performed with a common A'-factor for both Glasma and BFKL graphs A'ciasma = 
A'bfkl = 1.5. The bands in the plot correspond to different choices in Qq for the different 
track selections which are constrained to reproduce the A^"^™*^ in these windows. What 
these are can be deduced from the x-axis of Fig. ([2]) as discussed previously. We see that 
the agreement of the theory curves to the data in the higher multiplicity windows is quite 
good. 

However, there is a significant discrepancy in the two lowest multiplicity windows, es- 
pecially at higher pt- The reasons are two fold. Firstly, our formalism for both Glasma 
and back-to-back contributions is less valid for these multiplicities, where more peripheral 
impact parameters in the proton are accessed. However, a more important reason may be 
that the experiments require at least two particles in each bin to extract a di-hadron signal. 
Our theory computation (in absence of Monte-Carlo simulations that are challenging for in- 
terference graphs) imposes no such restriction. In lower multiplicity windows, and at higher 
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Pt, this is likely to provide a significant correction, just as seen in the theory comparison to 
data. 

Results for the coUimated yield versus A0 matrix in p'rP^^p^^ for the high multiplicity 
-^trk™*' — window are shown in Fig. (|4]). The difference to the plot shown previously in 
[l| is the common K factor for both BFKL and Glasma graphs of i^" = 1.5. Further, the 
Qo's corresponding to the high multiplicity window are higher. As noted earlier, this change 
was driven by a better understanding of the normalization of the different experiments, in 
this case CMS, subsequent to our previous papers. The QqS chosen are seen in Fig. ([2]) 
to better represent the centrality classes in multiplicity for the detailed matrix comparison. 
Given these noted changes, the agreement with the data over nearly 400 data points is 
remarkably good. While the Glasma signal is small in the data, the back-to-back correlation 
is significant, and it is striking that the BFKL graph captures its systematics so well 40 . 
As we argued in [l[, the multi-regge (QMRK) 2 — )■ 4 di-jet contribution without BFKL 
evolution between triggered gluons shows a significantly larger collimation, a conclusion 
that remains unchanged. The data in our view is demonstrating decorrelation of the back- 
to-back di-hadron signal due to QCD evolution a la BFKL between the triggered hadrons. 

We now turn to a discussion of a comparison of our Glasma-|-BFKL framework to first 
data from proton- lead collisions at the LHC. In [2|, we made a comparison to first data from 
the CMS collaboration [4]. In this paper, we will revisit that comparison for the reasons 
articulated above. Further, we will make quantitative comparisons with data from ALICE 
and ATLAS within their distinct experimental acceptance. As noted, a comparison of these 
is discussed in the appendix. Though first PHENIX data on very central deuteron-gold 
data has been presented at a conference [3], we are unable to make a direct comparison 
because the quantity presented, unlike the coUimated yield, is sensitive to the combinatorial 
background. We do however make a prediction for the coUimated yield. 

To simulate the p+Pb collision, we vary Qq at the initial rapidity scale in the proton and 
lead nuclei. All other parameters are the same, with the exception of -ft'ciasma = -^bfkl = 1, 
as opposed to 1.5 for proton-proton collisions. There is no reason a priori that K factors 
should be the same in the p+p and p+Pb case. It is conceivable that the dense-dense 
factorization used here is more applicable in the latter with smaller higher order corrections. 
In our treatment, the proton is varied in multiples AT^art""^ of the "minimum bias" value 
of Ql = 0.168 GeV^ to simulate proton- lead collisions that select more central impact 
parameters in the proton, where the gluon density is considerably higher than the gluon 
density for the median impact parameter corresponding to minimum-bias events. On the 
lead side, as noted previously, the initial saturation scale in lead is Ql = N^^^.^ ■ 0.168 GeV^, 
where A^p^^ denotes the number of participants (color charge probed) in the lead nucleus. 

SpecificaUy, to compare to the nearside coUimated yield data in Fig. (|5]), the lower and 
upper curves correspond to Ql values that are respectively five to six times the minimum 
bias value. As seen from the x-axis of Fig. ([2]), they provide a reasonable estimate of the 
^trk™°'s one can estimate contributing to the centrality cut A^°r^"^'^ > 110. By the same 
logic, for the same centrality cut in proton-lead coUisions, the guidance afforded by Fig. ([2]) 
suggests {Npl°l°'^,N^^^^) =(3,22) (upper curve) and (4,14) (lower curve). These are of course 
estimates, but we have checked that small variations of these do not significantly widen the 
uncertainty band. 

We see from Fig. that the nearside yield from the Glasma graphs, within theoretical 
uncertainties, is able to account for the pr = p%^^ = Pt^ dependence of the CMS mea- 
surement for the high multiplicity window in both proton-proton and proton-lead collisions. 
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Associated Yield 




FIG. 5. The px [p^^ = p^^) dependence of the associated yield in proton-lead and proton-proton 
collisions. The data here are for N°^^^'^ > 110. With the guidance from Fig. ([2]), the proton-lead 
centrality band correspond to (-^part°°,-^]^jjjt) of (3,22) (upper curve) and (4,14) (lower curve). The 
proton-proton curves correspond to N^l"^"^ = 5 and 6. See text for further explanation. 



In particular, it naturally explains the factor of six enhancement in proton-lead to proton- 
proton collisions in the pt = 1 — 2 GeV window where the yields are the largest. While some 
final state rescattering cannot be ruled out, we believe it will be difficult to account for this 



large factor in hydrodynamic models [41[. In our picture, the signal is due to the quantum 
interference of the unintegrated gluon distributions in the projectile and target. As shown 
in [2I, the contribution of this overlap scales approximately as A^t°k^'^'^ ' ^part; a result of 
which one obtains the trajectories shown in Fig. ([2]). 

A more detailed comparison of the CGC EFT framework to data is obtained in Fig. (Q. 
Here the CMS proton-lead data for the collimated yield as a function of A0 are shown for 
a number of windows in A^^°f'"'^ and in pt = p^'^ ~ Pt"^- The saturation scales for each 
N"^^'^'^ centrality window are estimated with the guidance from Fig. As in the proton- 
proton case, we see that the agreement is quite good, especially in the higher multiplicity 
windows. Again, as in the p+p case, we see the most significant deviation from data is 
the underprediction in the lower multiplicity bins at higher pt- We believe this to have the 
same underlying cause as in the proton-proton case; a trigger-bias in the experiment where 
only events containing at least two hadrons in the pt windows of interest are included in 
the averaging. At lower pt and higher multiplicities, where the yield of charged particles is 
larger, this effect becomes insignificant. 

Recently, subsequent to our paper j2|, papers on long ran ge r apidity correlations in the 
LHC proton-lead data at = 5.02 GeV from both ALICE [l^ and ATLAS [l8| collabo- 
rations have appeared. The former took data on very central events in the peudorapidity 
window \Ari\ < 1.8, while the later has a pseudorapidity acceptance 2 < |A?7| < 5, closer to 
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FIG. 6. Long-range (2 < |Ar/| < 4) per-trigger-yields {l / Nin^cP N / d/S.(j)^ of charged hadrons as a 
function of |Ai;^| in different px and multiplicity bins for pPb collisions at ^/s = 5.02 TeV with 
the CMS experiment. The pT selection corresponds to both particles in the pair. The centrality 
dependence of the theory curves are controlled by the choice of initial saturation scale in the target 
and projectile. These are i) N^^'""" < 35: {N^lT"" , ^p^rt) = (1> 3) (lower curve), (2, 6) (upper 
curve), ii) 35 < A^t°^'°^ < 90: (2, 6) (lower), (2, 12) (upper), iii) 90 < iVt°rk''''' < (2,14) (lower), 
(2,22) (upper), iv) 110 < N°f'''' < 150: (3, 22) (lower), (4, 14) (upper), and v) iV°f''°'= > 150: (4, 
16) (lower), (4,20) (upper). 



the CMS acceptance. (For a detailed discussion of the different acceptances and normahza- 
tion procedures of the different experiments, we refer the reader to the appendix.) A novel 
feature first introduced by ALICE, and shortly thereafter by ATLAS, is the subtraction 
of the associated yield per trigger in peripheral proton-nucleus from the same quantity in 
more central collisions. In our framework of Glasma graphs + BFKL graphs independently 
contributing to the yield per trigger, this procedure is especially valuable because it is equiv- 
alent to isolating the Glasma graph component j42|. As we checked previously, the BFKL 
di-jet per trigger contribution is very weakly centrality dependent [l|. 

In Fig. d?]), we show a comparison of results in our framework to the ATLAS data in sym- 
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FIG. 7. Left: Data from the ATLAS collaboration for the associated yield versus lS.(j) for the central 



and peripheral events identified in 13] ■ The blue theory curve corresponds to the comparison of 
Eq. ([5]) for the min-bias Q'^{proton) = 0.168 GeV^ and N^^^^^ = 3 to the peripheral collision data. 
The other theory curve is a comparison of Eqs. ([I|)+ ^ to the central data. The band corresponds 
to initial saturation scales, in our notation, of (4,14) (lower curve) and (3,22) (upper curve). Right: 
direct comparison of the central minus peripheral subtraction performed by ATLAS to the Glasma 
graph ([T]) contribution. 



metric px windows. The two data sets at higher values are also available from the CMS 
experiment; as noted, ATLAS has a larger Ar/ acceptance relative to CMS. In the left plot, 
we show a comparison of the ATLAS peripheral data to the BFKL contribution+ Glasma 
contribution. The Glasma contribution is negligible in peripheral events. We see that 
^part""^ = 1 on the proton side and N^^^^ = 3 on the lead side give a good description of the 
data. This adds confidence to our fits to the CMS data since the quoted values for Qq are 
indeed of the order of what one would expect as typical values in peripheral collisions. The 
central ATLAS data should be compared to the sum of Glasma + BFKL contributions, and 
we see again we get good fits for the band (3,22) (upper curve), (4,14) (lower curve), as in 
the CMS comparison. The agreement is quite reasonable though the theory curves slightly 
under-predict the data in the plots of the net associated yield on the nearside (left plots 
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FIG. 8. Per-trigger-yield I / Ntrigd^ N / dA(j) versus |A(/>| for asymmetric p^'^, p^^. Red squares 
are the ATLAS data for central p+Pb events, blue circles are for peripheral p+Pb events. The 
blue curves show the results of the BFKL contribution for (A^p^°*°", A''^jjj.^)=(l,3). The gray band, 
representative of central collisions, corresponds to (4,14)-lower curve and (3,22)-upper curve. 
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FIG. 9. Comparison of Glasma graphs, Eq. ([TJ, to ALICE data for the cohimated yield per trigger 
per Ar] for central (0-20%) events with (60-100%) peripheral events subtracted versus Ac/). The 
band is from a computation with (3Q(proton) = 0.336 GeV^ on N^^^.^ = 12 (lower) and Np^j.^ = 14 
(upper). 

of Fig. ([7])). We attribute this to an artificial anti-collimation of the BFKL contribution 
at small A0. After the jet subtraction performed by ATLAS the agreement between the 
Glasma Graphs and ridge yield is restored as demonstrated in the right plots of Fig. ([7]), 
where the "di-jet" subtracted yield is compared to our Glasma graph. The overall agree- 
ment is quite reasonable suggesting a consistent interpretation of the power counting in our 
"initial-state" framework for multi-particle production and the experimental observations. 

The ATLAS collaboration has presented data for proton-lead collisions in a wide range 
of Pt^^Pt'^ windows for central and peripheral collisions. In addition to the symmetric 
Pt = P^T^ = Pt^ windows have already been presented by CMS for a similar acceptance, we 
consider in Fig. ([8]) a comparison to asymmetric p^'^ ^ p^'^ windows. We have also included 
in our comparison the data for p'^'^, where the lower range of the window is pt = 0.5 
GeV, which one might consider quite low p^'s for our framework. For the asymmetric 
windows shown in Fig. ([H]), the agreement on the awayside is quite good for both central 
and peripheral events, though arguably underpredicting the central events. However, on 
the nearside, the Glasma contribution is significantly lower than the data for central events, 
especially for the larger p^^^^ windows. Of all the p^^^^^p^^ windows presented by ATLAS, 
the asymmetric 3 < p^'^ < 4 GeV windows provide the worst comparison to the nearside 
Glasma computation. On the theory side, we have not attempted any fine tuning with small 
systematic adjustments of the K factors and Ql values. At these higher pt values the same 
trigger bias, as discussed earlier, has to be taken into consideration. It will be interesting 
to see if this discrepancy persists with the additional data that is anticipated to be released 
soon, and if this discrepancy in asymmetric windows is also seen by CMS and ALICE. 

Data on proton-lead collisions from the ALICE collaboration is shown in Figs. ^ and 
(ITU]) . The acceptance of the ALICE experiment is distinctly different from CMS or ATLAS, 
covering IAt^I < 1.8. For small [Ar^l < 1, there is a nearside A0 di-hadron correlation 
that can be attributed to jet fragmentation. If this short range component can be safely 
subtracted, one can look for a contribution that is long range in rapidity. Further subtraction 
of the awayside "di-jet" contribution, whose yield per trigger is observed to be weakly 
centrality dependent, will, as discussed, reveal the Glasma graph contribution. The Glasma 
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FIG. 10. Comparison of the (60-100% centrality subtracted) data for both nearside and awayside 
integrated colhmated yields from different centrahty classes for "off-diagonal" windows (2 < p^'^ < 
4 GeV: i) 0.5 < pf?^ < 1 GeV, ii) 1 < < 2 GeV, and 1 < p^'^ < 2 GeV, 0.5 < pf?^ < 1 Gev) 
and symmetric windows (p^^^ = p'^'^ = 0.5 — 1,1 — 2,2 — 4 GeV). The band for the 0-20% centrality 
class is the result for (2,12) (lower curve) and (2,14) (upper curve), 20-40% centrally: (2,4) (lower), 
(2,6) (upper) and 40-60% centrality: (1,3) (lower), (1,4) (upper). 



graph contribution is nearly rapidity independent (for |A?7| ~ l/as 5 units of A?]), so 
one expects such a contribution to be present in the ALICE acceptance. In Fig. ([9]), we 
show the comparison of the Glasma graph computation to the central (0 — 20%) minus 
peripheral (60 — 100%) yield from ALICE. The agreement is quite good. The overshoot on 
the awayside is sensitive to smearing of the back-to-back contribution and possible systematic 
uncertainties due to the peripheral jet subtraction. 

We now turn to the description of the integrated nearside and awayside ALICE data in 
Fig. (fTOl) . We see that for nearly all the symmetric and asymmetric, near and awayside pt 
windows, and all centrality classes listed in the caption of Fig. ffTOl) . there is good agreement 
with the ALICE data. There is a slight overshoot of the data in the lowest pt window for the 
most central collisions and for one of the windows in the 40-60% centrality range. The ALICE 
paper 17|] also quotes values for the V2 and flow moments. In our framework, extraction 
of these moments also depends on the combinatorial Acf) independent background, which 
varies from one pt window to the next. Since there are many possible QCD contributions 



to the combinatorial background [43|-|46|, an estimate of this quantity is less reliable in 
our framework, though in principle feasible in future. For the connected graphs we have 
considered so far, there can be no ^3 contribution because these Glasma contributions are 
symmetric about A0 = 7r/2. It remains an open question whether the full set of connected 
graphs can produce a small collimation due to weak final state effects. In any event, 
a large component would be challenging for our framework. Our take on the ALICE 



f3 data presented in [17|| is that the effect observed is sensitive to systematics of the jet 



subtraction-see also the comment in footnote IIVI Further data from the 2013 p+Pb run 
should help clarify these issues significantly. 

Last but not least, in Fig. (fTTi) . we show predictions for the collimated yield from a 
correlated di-hadron signal that is long range in Ar/ for RHIC energies. It is clear from 
Fig. ([2]) that a collimated signal is feasible at the significantly lower energies of deuteron- 
gold collisions at RHIC as long as very central events are triggered on. The predictions for 
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FIG. 11. Predictions for the associated yield versus A(p (left figure) and the A(p integrated associ- 
ated yield p^'^ (right) figure from Glasma graphs for RHIC energies. The result from Fig. [10] for 
the ALICE experiment is show for comparison. 

different centrality classes and the functional dependence of the coUimated yield on p^^^ are 
shown in Fig. f|TT]) -the ALICE prediction is shown in comparision. A cautionary aspect of 
comparisons to RHIC energies is that the higher pt windows and Ar] > correlations are 
sensitive to a; > 0.01 values in the nucleon and nuclear wavef unctions. The CGC framework 
should break down at these large x values. For modest plp^,p^^ windows at RHIC, where the 
physics is still weak coupling and x is still small, one expects important information due to 
the widely different energy to help constrain our picture of long range rapidity correlations. 
Very recently, first data from these correlations in deuteron-gold collisions at y/s = 200 GeV 
from the PHENIX experiment were presented [19] . The analysis was very similar to that of 
the ATLAS experiment. A significant V2 is observed, while is consistent with zero. As 
we noted previous, due to the sensitivity of the result to the combinatorial background, a 
direct comparison of results similar to Fig. (ITTj) to the data cannot yet be achieved. We 
expect that results for the coUimated associated yield per trigger will become available soon 
facilitating a more direct comparison to our predictions for RHIC energies. 



IV. SUMMARY AND OUTLOOK 

In this third paper in the series, preceded by [H, 0], we developed significantly the com- 
parison of the CGC EFT to large At] di-hadron correlation data from the LHC experiments. 
The different normalizations of data taken by the different experiments were taken into ac- 
count and a consistent treatment of CMS, ALICE and ATLAS data was presented - the last 
two being discussed for the first time. In the CMS case, we presented a reanalysis relative 
to [H, 0| both for the proton-proton and proton-nucleus data. In the former case, we showed 
for the first time detailed comparisons of theory to data from a number of N^^^^'^ windows. 

The agreement of theory with data in proton-proton and proton-lead collisions over a very 
wide range of p^rP^,p^^ windows, centrality class and Ar] acceptance, is quite spectacular 
and lends strong support that a) gluon saturation is being seen in these experiments, and 
b) the data are sensitive to systematics of renormalization group evolution of unintegrated 
gluon distributions (UGD) in the description of Glasma and "back-to-back" graphs. It 
is remarkable that gluons widely separated in rapidity show a A0 ^ coUimation that 
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is sensitive to detailed dynamical features of the theory. Because the Glasma graphs are 
quantum interference graphs, which have different structures in the amplitude and complex 
conjugate amplitude, the result is a particular form of gluon entanglement, unique to QCD. 
If this picture is confirmed by further data, it provides an impetus for further developments 
in theory to better understand the properties of saturated gluon states in the nucleon and 
nuclear wavef unctions, providing an important window to hadron structure and dynamics 
at high energies. 

An alternative scenario for the dynamics of long range rapidity correlations is from flow 
resulting from strong final state rescattering jiH . We beheve this possibility is 



strongly disfavored on both conceptual and phenomenological grounds-these issues will be 
addressed elsewhere. Our perspective based on the discussion in this paper is that the onus 
is now on models based on final state scenarios to present a quantitative description of data 
that is as transparent and competitive to the same degree as the one presented here-for 
instance, the matrix of data of the coUimated yield versus A0 for the various experiments. 

Further data from the LHC will no doubt provide clarity and challenge either or both 
scenarios. Regardless, unless the data presented thus far changes significantly, the results 
of our detailed analysis suggest that initial state effects must play an important role even if 
other QCD effects come into consideration in the description of these striking experimental 
phenomena. 
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APPENDIX: EXPERIMENTAL COVERAGE AND RELATIVE NORMALIZA- 
TIONS OF DI-HADRON DATA IN THE LHC EXPERIMENTS 

In this appendix, we discuss the relative experimental acceptances and normalizations of 
the different LHC experiments. The rapidity acceptance of these are 
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The effect of different normalizations between the experiments can be gauged by consid- 
ering single inclusive and double inclusive charged hadrons obtained from Eqs. (|9]) and (fTOj) . 
In this case, 



and 

Jvmin ^Amin J P^'^ ^ J q^'^ ^ ^{Vp^Vq) dripdrigd^pxd'^ 



(13) 
where 

A {r]p, r]g, (t)p, (pq) = 5{(t)p - (pq - A(l))9 {\r]p - r]g\ - Ar/min) (Ar/max - \Vp - VqD ■ (14) 

In Eq. f lT3|) . the ATLAS experiment does not weight their signal by the background nor 
do they calculate their result per unit Arj. For ATLAS, we have therefore B = 1. For CMS 
and ALICE, 

\Vp - Vq\ 



BcMS = -BalICE = 2|A?7max - At^i, 



I ^max ^min | 



(15) 



To get a sense of what these different normalizations entail, let us compute 
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One will then have 
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(18) 



There are additional subtleties: 



The CMS experiment presents the per-trigger-yield differential in A0 between to 
TT. There is a corresponding yield from to — vr which is not shown nor included in 
the contribution of the associated yield (AY). Both ALICE and ATLAS include the 
contribution from to — vr when computing the AY thus making their results a factor 
of two larger. When the ATLAS plots are shown from to vr they fold over the results 
from —IT to making the differential yield a factor of two larger. ALICE plots the 
differential yields over the full phase space of — tt to vr therefore making a correction 
unnecessary. 
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ALICE has a cut on their integrals such that p^'^ < pf?'^ for symmetric windows 
the signal will be a factor of two smaller than for asymmetric windows. 



With these considerations, we obtain 
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